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Abstract In the present paper we give a historical account -ranging from classical 
to modern results- of the problem of rolling two Riemannian manifolds one on the 
other, with the restrictions that they cannot instantaneously slip or spin one with 
respect to the other. On the way we show how this problem has profited from the 
development of intrinsic Riemannian geometry, from geometric control theory and 
sub-Riemannian geometry. We also mention how other areas -such as robotics and 
interpolation theory- have employed the rolling model. 



1 Introduction 

Differential geometry has been inextricably related to classical mechanics, since 
its very conception in the 18th century. As a matter of fact, back in the days, this 
area of research was referred to as rational mechanics. The basic idea of this point 
of view is reasonably simple: to a given mechanical system M, one can associate 
a differentiable manifold M in such a way that each possible state of the system 
corresponds to a unique point in M. In this way, each possible velocity vector of M 
at a given configuration is represented as a tangent vector to M at the corresponding 
point. The classical dictionary goes as follows: 
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1 . Physical data (such as masses, lengths, etc.) of elements in M induce a Rieman- 
nian metric in M representing the kinetic energy. 

2. Linear restrictions imposed on the positions of M (or that can be integrated to 
such) translate to submanifolds of M. 

In the late 19th century, physicists noted there were plenty of mechanical systems 
not considered by the above dictionary. These systems were named non-holonomic, 
opposed to holonomic systems which are defined in the second point of the dictio- 
nary above. A mechanical system M is non-holonomic if its dynamics has linear 
restrictions that cannot be integrated to constraints of the position. For various ex- 
amples and a brief historical bibliography, we refer the interested reader to the sur- 
vey JH). A well-known early example of these systems is the sphere rolling on the 
plane without sliding or spinning, studied (with some variants) by S. A. Chaplygin 
in the seminal works |[T6l [171 . Our aim in this paper is to give a general look at 
some of the most important breakthroughs in mathematics that gave us some un- 
derstanding of the generalized version of this system consisting on two Riemannian 
manifolds M and M of the same dimension rolling one against the other, not al- 
lowing instantaneous spins or slips. Nowadays these systems are often studied in 
connection to sub-Riemannian and Riemannian geometry fl3l l48l and geometric 
control theory J5] . 

The structure of the paper is the following. In Section [2] we recall two major 
players in the study of the mechanical system described above and early differen- 
tial geometry: S. A. Chaplygin and E. Cartan. Chaplygin studies for the first time 
the problem from a mechanical point of view and finds first integrals of motion in 
different situations. Cartan's development and his celebrated "five variables" paper 
were not evidently connected to the rolling model at the time of their publication, 
see iflOl . nevertheless we present them from our point of view. In Section [3] we 
briefly present Nomizu's breakthrough introduction of the dynamics of rolling in 
higher dimensions, through embedded submanifolds of Euclidean space and its rela- 
tion to Cartan's development. In Section|4]we present how the problem was brought 
back to life when control theory sees in differential geometry a useful tool to treat 
the controllability issue of the rolling model in two dimensions and some geometric 
consequences of optimality conditions. Section [5] surveys how the higher dimen- 
sional rolling system was re-discovered and how it appears naturally in geometric 
interpolation. Finally in Section [6] we present the latest results that have been ob- 
tained concerning the controllability of the system and its symmetries. We conclude 
with a brief discussion on some generalizations and open problems. 



2 The early years: Mechanics and the new differential geometry 

The first time the problem of a ball rolling on the plane was considered as worthy 
of study was in the seminal papers of S. A. Chaplygin lfl6l[T7l . one of the fathers 
of non-holonomic mechanics. The results were considered surprisingly difficult at 
the time, and for [ 16 1 Chaplygin won the Gold Medal of Russian Academy of Sci- 
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ences. The main results he obtained were first integrals of motion for the system in 
several geometric situations. Even these seemingly elementary problems contains 
unexpected difficulties and bottlenecks when trying to obtain closed formulae for 
the dynamics. As stated in |[T6ll . after observing that the differential equation of the 
dynamics can be integrated in quadratures. Essentially at the same time, E. Cartan 
was developing his coordinate-free differential geometry. With this new language 
he was able to propose and study many problems, most often related to the search 
of invariants of geometric systems. In this survey, we will only focus in two of his 
many ideas: the search for invariants and symmetries for control systems with two 
controls and five degrees of freedom, and the definition of affine Riemannian holon- 
omy through the development of a curve. Both of this ideas will appear several other 
times in this survey. 



2.1 Chaply gin's ball 

In the year 1897 the work lfl6] written by S. A. Chaplygin was published. This 
papers is one of a series of research articles in which Chaplygin analyzed non- 
holonomic systems. Also of particular relevance to this survey is another paper ifTTl . 
In particular he was interested in studying first integrals and equations of motion for 
different systems of rolling balls. 

To illustrate his results, Chaplygin was able to find an integral of motion for the 
system of a homogeneous small ball of mass m\ and a homogeneous sphere of mass 
ni2, in which the ball rolls without slipping inside the sphere. We will think of the 
dynamics occurring in Euclidean 3-space. Let O be the center of the sphere, let G be 
the center of the moving ball and A the point of contact between the two. Introducing 
the quantities a = dist((9, G) and b = dist(<9, A), then one has the integrals of motion: 

2 / dxi dy { \ t Jb .\( da dp\ 

Y m yt — - -x,— +M --1(8 a— = const. 

g V dt dt J \a )\ dt dt ) 

Where A = (a, 6, 7) with respect to a fixed frame OX'Y'Z', and the points G = 
(jci,yi,zi) and O = (x2,y2,Zi) with respect to a moving frame AXYZ, with axes at 
all times parallel to those in OX'Y'Z'. Additionally M = nt\ +1112 denotes the mass 
of the system. 

The equations of motion are complicated and it serves little purpose to write them 
down here. Nevertheless, there is an interesting historical remark at this point. After 
arriving at a very complicated differential equation to describe the dynamics of the 
system, Chaplygin observes it can be written in the form 



- 7 + v<£(C)+f(C) -0, 
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for some appropriate functions <P and f after a series of changes of variables. He 
then ventures to say 

[. . . ] and, therefore, can be integrated in quadratures. We will not write out these quadratures 
since they are rather cumbersome. 

As far as we know, the integration of differential equations connected to the problem 
of rolling balls is still an area of active research, see for example [ 13 1. 



2.2 Cartan's "five variables" paper 

A rank / vector distribution D on an n-dimensional manifold M or (Z,«)-distribution 
(where / < n) is, by definition, an /-dimensional subbundle of the tangent bundle 
TM, i.e., a smooth assignment q i-> D\ q defined on M where D\ q is an Z-dimensional 
subspace of the tangent space T q M. Two vector distributions D\ and D2 are said to be 
equivalent, if there exists a diffeomorphism F : M — » M such that F*Di \ q = D^ptq) 
for every q £ M. Local equivalence of two distributions is defined analogously. 

Cartan's equivalence problem consists in constructing invariants of distributions 
with respect to the equivalence relation defined above. A seminal contribution by E. 
Cartan in [14] was the introduction of the "reduction-prolongation" procedure for 
building invariants and the characterization for (2, 5) -distributions via a functional 
invariant (Cartan's tensor) which vanishes precisely when the distribution is flat, 
that is, when it is locally equivalent to the (unique) graded nilpotent Lie algebra f) 
of step 3 with growth vector (2,3,5). 

In the same paper, Cartan also proved that in this system there is hidden a real- 
ization of the 14-dimensional exceptional Lie algebra 02- To explain where does it 
appear, let us recall that an infinitesimal symmetry of an (l,n)— distribution D is a 
vector field X E VF(M) such that [X,D] C D. Now consider the (unique) connected 
and simply connected nilpotent Lie group H with Lie algebra h. The two dimen- 
sional subspace of f) that Lie generates it, can be seen as a (2,5)— distribution on H. 
In general, a (2,5)— distribution that is bracket generating is nowadays known as a 
Cartan distribution. In this setting, the following theorem takes place. 

Theorem 1 (Cartan 1910). The Lie algebra of symmetries of the fiat Cartan dis- 
tribution is precisely Q2> an d this situation is maximal, that is, for general Cartan 
distributions the dimension of the Lie algebra of symmetries is < 14. 

Moreover, Cartan gave a geometric description of the flat G2-structure as the 
differential system that describes space curves of constant torsion 2 or 1/2 in the 
standard unit 3-sphere (see Section 53 in Paragraph XI in lfi"4"l .) 

The connection between this studies by Cartan and the rolling problem comes 
from the fact that the flat situation described above occurs in the problem of two 
2-dimensional spheres rolling one against the other without slipping or spinning, 
assuming that the ratio of their radii is 1 : 3, see |[T2| for some historical notes and a 
thorough attempt of an explanation for this ratio. In fact, whenever the ratio of their 



The rolling problem: overview and challenges 



5 



radii is different from 1 : 3, the Lie algebra of symmetries becomes so(3) x so(3), 
thus dropping its dimension to 6. A complete answer to this strange phenomenon 
as well as a geometric reason for Cartan's tensor was finally given in two remark- 
able papers ll52l l53l (cf. also [4|), where a geometric method for construction of 
functional invariants of generic germs of (2,«)-distribution for arbitrary n > 5 is de- 
veloped. It has been recently observed in [5 1 that the Lie algebra of symmetries of a 
system of rolling surfaces can be 02 in the case of non-constant Gaussian curvature. 



2.3 Cartan 's development 

E. Cartan in |fT31 defined a geometric operation, that he called development of a 
manifold onto a tangent space, in order to define holonomy in terms of "Euclidean 
displacements", i.e., elements of E(«). In his own words: 

Quand on developpe l'espace de Riemann sur l'espace euclidien tangent en A le long 
d'un cycle partant de A et y revenant, cet espace euclidien subit un deplacement et tous 
les deplacements correspondant aux differents cycles possibles forment un groupe, appele 
groupe d'holonomie. 

An interpretation of this quote in terms of manifolds rolling follows naturally. For 
a given loop y: [0, t] — » M on an n dimensional Riemannian manifold M, one can 
roll M against the Euclidean space M" obtaining a new curve f: [0, t] — > K", called 
the development of y. By parallel transporting along y any orthonormal frame of 
T\y(mM, we obtain a rotation Ry £ 0(n). The fact that y is not necessarily a loop 
induces a translation Ty corresponding to the vector y(r) — y(0). We conclude that 
we can associate to y an element (Ry, Ty) of the Euclidean group of motions E(n). 
The subgroup Hol"^ (M) of E(n) consisting of all such (Ry, Ty) obtained by rolling 
along all absolutely continuous loops y is known as the affine holonomy group of 
M and the orthogonal part Hol(M) C 0(n) of it is the holonomy group of M. 

It is known that if M is complete and with irreducible Riemannian holonomy 
group, the affine holonomy group contains all translations of T \ X M, see 11371 Corol- 
lary 7.4, Chapter IV]. In other words, under the irreducibility hypothesis, the rota- 
tional part of the affine holonomy permits to recover the translational part, and this 
consists of all the possible translations in T\ X M. 

Perhaps something that might have been not expected by Cartan is that this con- 
cept of development would play a fundamental role in the definition of Brownian 
motion on a manifold, and the subsequent explosion of interest that stochastic anal- 
ysis in Riemannian manifolds has had in later decades, see (29]. For a long time, 
mathematicians have had the intuition that by rolling an n-dimensional manifold 
M along a given curve y(t ) in R" with the Euclidean structure, one would obtain a 
curve in M which resembles the original curve y(t), see ll27l . The main outstanding 
idea (as far as we know due to Malliavin) was to use Cartan's development through 
the orthonormal frame bundle and Wiener's measure, see ll50l . 

The idea of how to define Brownian trajectories on manifolds is similar to the 
interpretation given above. Intuitively, one can draw a Brownian path B(t) in M", 
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and then one can consider the system of M rolling against R" following the path 
B(t). The precise definition uses a less regular version of Cartan's development and 
parallel transport. 

This naive notion allows one to recover the Laplace-Beltrami operator Am of 
the manifold. It is often interpreted as if Brownian paths are the "integral curves" 
for Am- Of course this assertion lacks of mathematical precision, but it introduces 
the idea that second order differential operators induce "diffusions" on the mani- 
fold. This point of view has been exploited significantly in the study of stochastic 
differential equations on manifolds, see Q. 



3 A "forgotten" breakthrough 

An important contribution to the understanding of the problem of rolling without 
slips or spins came to light in the paper [45 1 by K. Nomizu. His aim was to give a 
mechanical interpretation of certain differential geometric invariants using this sys- 
tem. He mainly focuses in submanifolds of M. N with the usual Euclidean structure, 
and so will we along this section. 

He begins with a simple general consideration: as a motion occurring in a Eu- 
clidean space M. N without deforming objects, a rolling can be seen as a curve in the 
Euclidean group E(N), that is a function [0, t] 9 f H> /> € E(JV) given by 

*-(o?)- (1) 

where / = Id is the identity matrix of (N +1) x (N + 1), Q € 0(N) and c, 6 R^. 
He calls such types of curves 1 -parametric motions. 

For a given 1 -parametric motion {ft}, he observed that there is a natural time- 
dependent vector field X t associated to it. For an arbitrary point y G R we define 

{X t ) y :— , where x — f~ l (y)- Using equation (Q~|), one can see that (X t ) y — 

u u=t 

S,y + v t , where S, = ^C^ 1 G o(N) and v, = -S t c, + ^ € R N . The corresponding 
element of the Lie algebra e(N) 

df r i = (S,v,\ (2) 

dt h ^ooj {Z> 

is called the instantaneous motion. Slips and spins can now be encoded in terms of 
the vector field X t and the instantaneous motion. 

Definition 1. The instantaneous motion (f2]i is called an instantaneous: 

• standstill if S t = and v t = 0, 

• translation if S t = and v t ^ 0, 

• rotation if there exists a point yo € M. N such that (X t ) yo = and St ^ 0. 
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With this at hand, it is possible to define rolling without slipping (skidding in 
Nomizu's terminology) nor spinning between M",M" «-> M. N . 

Definition 2. Let {/,} be a 1 -parametric motion such that f(M) is tangent to M at 
a point y t E M. Assume that (X t ) yt = and S t ^ 0. The motion / ( is a rolling if for 
any pair of tangent vectors X,Y E T yt N 

(S t (X),Y}=0, (3) 

and for any pair of normal vectors U, V € Ts~M 

(S,(U),V)=0. (4) 

An equivalent way of stating conditions (O and is that S t maps T yt M to T^M 
and also maps Tr-M to T yt M, 

This definition allowed Nomizu to find a very concrete realization of Cartan's 
development. For the case of surfaces rolling on the plane, his result reads 

Theorem 2 (Nomizu 1978). Let x t be a smooth curve on a surface M which does 
not go through a flat point of M. There exists a unique rolling {f} of M on the 
tangent plane E at xo such that y t — ft(x t ) is the locus of points of contact. The 
curve yt is the development of the curve x t into E. 

As a consequence of this result, Nomizu noticed that there is a natural kinematic 
interpretation of the Levi-Civita connection for a surface M, coming from the rolling 
formulation: a vector field U(t) along the curve x t is parallel with respect to the 
Levi-Civita connection of M if and only if C t (U(t)) is a constant vector for all t. 

As a matter of fact, he was able to extend this result to higher dimensions and 
gave conditions under which rollings exist in terms of the shapes of the submani- 
folds, that is, in terms of both intrinsic and extrinsic data. 

For reasons unknown to us, this paper seems to have been forgotten over the 
years. Nomizu's definition of higher dimensional rolling is equivalent to Sharpe's 
one in Subsection 15. H and many of his observations have been rediscovered in [48 
Appendix B]. Nevertheless, there is no reference to the paper P31 in Sharpe's book. 



4 Revival: The two dimensional case and robotics 

The aim of this section is to put in context the study of the rolling model for the 
case of two dimensional manifolds, and how they appeared naturally in problems of 
sub-Riemannian geometry, robotics and geometric control theory. 
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4.1 Rigidity of integral curves in Cartan's distribution 

In the celebrated paper 0111 . R. Bryant and L. Hsu studied curves on a manifold 
Q of dimension n > 3 tangent to a (2, n)— distribution D. The idea was to analyze 
the space QD(p,q) of differentiable curves in Q connecting two points p,q £ Q and 
being tangent to D (called D-curves by them). The space Qo(p,q) is endowed with 
its natural C 1 topology. The idea that D-curves can be "rigid" plays a fundamental 
role in their paper. 

Definition 3. A D-curve y: [0, t] — > Q is rigid if there is a C 1 -neighborhood % of y 
in f2 D (y(0), y(f)) so that every j\ £ is a reparametrization of y. We say that y is 
locally rigid if every point of / = [0, t] lies in a subinterval J C / so that y restricted 
to J is rigid. 

Their main result goes as follows. 

Theorem 3 (Bryant & Hsu 1993). Let D be a non-integrable rank 2 distribution 
on a manifold Q of dimension (2 + s) > 3. Suppose further that the distribution 
D\ = [D,D\ (which has rank 3) is nowhere integrable. Then there always exist D- 
curves that are locally rigid. 

They give a more precise description of such curves in terms of projections of char- 
acteristic curves in a dense subset of the annihilator of D\, but stating it precisely 
would not serve the purposes of this exposition. 

For us, the most relevant part of their work is their section on examples, in par- 
ticular their study of systems of Cartan type and of rolling surfaces. 

Recall that a bracket generating (2,5)— distribution is said to be of Cartan type. 
In other words D is a Cartan distribution if D\ has rank 3 and D2 = [D\ , D] has rank 
5. As a consequence of Theorem|3] they observe that there is exactly a 5-parameter 
family of locally rigid D-curves. In fact they briefly discuss a remarkable geometric 
behavior occurring in this situation: if M is connected, then any two points of M can 
be joined by a piecewise smooth D-curve, whose smooth segments are rigid. 

After all these observations, they devote themselves to the analysis of two 
oriented surfaces M and M endowed with Riemannian metrics rolling one over 
another without slipping or twisting. Let F and F be the oriented orthonormal 
frame bundles of M and M. Bryant and Hsu considered the "state space" mani- 
fold Q = (F x F)/SO(2), where SO(2) acts diagonally on the Cartesian product. 
An element in Q is a triple (x,x;A), where x £ M, x £ M and A : T X M — > T^M is an 
oriented isometry. Their formulation is as follows. Consider a curve y: [0, t] — > Q 
given by y(f) = (x(t),x(t);A(t)), then the no-slip condition reads A(t)(x(t)) —x(t). 
The no-twist condition requires some more care. Let ei,/i : [0, t] — > TM be a par- 
allel orthonormal frame along the curve x(t ) and let 

e 2 {t)=A{t){e l {t)), f 2 {t)=A{t){f l {t)), 

be the orthonormal frame along x(t) obtained via A. The rolling has no-twist when- 
ever the moving frame e2,/2 is also parallel (along x). 
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An important insight for the problem was expressing the no-twist and no-slip 
conditions in terms of a (2,5)— distribution D on Q. Let a,\ , (X 2 ,OC 2 i be the canonical 
1 -forms ofMonf and similarly )3i,j3 2 ,j3 2 i forM, see [49 1. Recall that these forms 
satisfy the so-called structure equations 

dai = a 2 iAa 2 , dj3i = foi Aj3 2 , 

da 2 = -a 2 iAai, dfh. = -j8 2 i Aj3i, 

da 2 \ = rca\ Aa 2 , dfoi = kpi Aj3 2 , 

where re and £ are the Gaussian curvatures of M and M respectively. With all of 
this, one can consider the distribution D on F x F defined by the Pfaffian equations 

ai - j8i = a 2 - j3 2 = a 2 i - j8 2 i = 0. 

The distribution they were looking for corresponds to the "push-down" image of D 
under the submersion F x F — > Q. A smooth curve y: [0, t] — > Q describes a rolling 
without slipping or twisting if and only if 7 is a D-curve. 

A remarkable fact is that the distribution D is of Cartan type whenever K — k ^ 0, 
which is an open set in Q. On this set, the corresponding 5 -parameter family of rigid 
curves describes the rolling of M against M following geodesies. 



4.2 Non-holonomy in robotics 

The traditional modeling of a mechanical system considers configurations (or states) 
of this mechanical system as points q of a smooth finite-dimensional manifold M, 
and the corresponding velocities q £ T q M are subject to locally independent con- 
straints in the Pfaffian form 

A(q)q = 0, (5) 

where A(-) is an m x n matrix of real-valued analytic functions, where m < n. Con- 
straints are said to be holonomic if their differential form given by (|5) is integrable. 
In this case, there exist integral submanifolds of dimension n — m that are invariant. 
If the constraints are not holonomic at some qoE M , then there will exist an integral 
submanifold containing q$ of dimension n — m + k with < k < m. The integer k is 
referred to as degree of non-holonomy. If k — m, the constraints, and by extension 
the system, are said to be maximally non-holonomic (see [44|). 

There is a more convenient way for control theory to describe the constrained 
system. If G(q) denotes a matrix whose columns form a basis for the annihilating 
distribution of A(q), then all admissible velocities q g A(q) ± C T q M can be written 
as linear combinations of the columns of G(q), 



n—m 

q = G{q)w = £ gi(q)wi, 



(6) 
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where w is a vector of quasivelocities taking values in IR" _m . When quasivelocities 
can be assigned values at will in time, functions can be regarded as control inputs 
of the driftless, linear-in-the-control, nonlinear system defined by ©. A physical 
actuator is associated to each control input , e.g. a motor for electromechanical sys- 
tems. The issue of non-holonomy of the original system, i.e. non-integrability of 
©, can be addressed by studying the distribution A spanned by the the vector fields 
gi's and more precisely the corresponding Lie algebra generated by the g/'s. If the 
system is maximally non-holonomic (or completely controllable), any two config- 
urations q and q' of its n-dimensional manifold can be connected along the flows 
of n — m vector fields. From an utilitarian engineers viewpoint, the latter definition 
may be rephrased as an «-dimensional non-holonomic system can be steered at will 
using less than actuators. This formulation underscores the appealing fact that de- 
vices with reduced hardware complexity can be used to perform nontrivial tasks, if 
non-holonomy is introduced on purpose, and cleverly exploited, in the device design 
(see 1441 "). 

Non-holonomy of rolling is particularly relevant to robotic manipulation, one of 
the main goals of which is to manipulate an object grasped by a robot end-effector so 
as to relocate and re-orient it arbitrarily, the so-called dexterity property. Dexterous 
robotic hands developed so far according to an anthropomorphic paradigm employ 
far too many joints and actuators (a minimum of nine) to be a viable industrial 
solution. Non-holonomy of rolling can be used to alleviate this limitation. In fact, 
while rolling between the surfaces of the manipulated object and that of fingers has 
been previously regarded as a complication to be neglected, or compensated for, 
some works (see, in particular, [1. 6. 18. 24. 39. 40] and the references therein) tried 
to exploit rolling for achieving dexterity with simpler mechanical hardware. 

Introducing non-holonomy on purpose in the design of robotic mechanisms can 
be regarded as a means of lifting complexity from hardware to the software and 
control level of design. In fact, planning and controlling non-holonomic systems 
is in general a considerably more difficult task than for holonomic systems. The 
very fact that there are fewer degrees-of-freedom available than there are configura- 
tions implies that standard motion planning techniques can not be directly adapted 
to non-holonomic systems. From the control viewpoint, non-holonomic systems are 
intrinsically nonlinear systems, in the sense that they are not exactly feedback lin- 
earizable, nor does their linear approximation retain the fundamental characteristics 
of the system, such as controllability (see $4\ ). 

The system of rolling bodies considered here differs substantially from the class 
of chained form systems or differentially flat systems (see Rouchon ll46l ). Con- 
sider, for example, the case such of the plate-ball system (i.e. a ball rolling on a 
plane without slipping or spinning), which is a classical problem in rational me- 
chanics, brought to the attention of the control community by Brockett and Dai |9|. 
Montana [42] derived a differential-geometric model of the rolling constraint be- 
tween general bodies, and discussed applications to robotic manipulation. Li and 
Canny 1 38 1 showed that the plate-ball system is controllable, and that the same holds 
for two rolling spheres, provided that their radii are different. 
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We close this subsection mentioning the beautiful works of Jurdjevic l34l [351 
who studied the problem of finding the path that minimizes the length of the curve 
traced out by the sphere on the fixed plane. It turns out that optimal paths also min- 
imize the integral of their geodesic curvature, so that solutions are those of Euler's 
elastica problem. For the higher dimensional cases of this problem, see Il36ll54l . 



4.3 Orbits and complete answer for controllability 

The point of view adopted by Bryant and Hsu was improved significantly by A. 
Agrachev and Y. Sachkov in |2| employing tools in geometric control theory. 

Two innocent, yet powerful, changes in perspective made the problem more ac- 
cessible for the application of the orbit theorem of Sussmann [51 1. These modifi- 
cations consist of rewriting the state space of the rolling and, most importantly, to 
prefer the use of vector fields (written in local coordinates) instead of differential 
forms (written without using coordinates). 

Let M and M be smooth two-dimensional connected oriented Riemannian sur- 
faces. The new version of the state space is given by 

Q = {A : T X M -> T$M\x £ M,x € M,A an oriented isometry}. 

It is an easy exercise to see that Q is indeed diffeomorphic to the manifold M intro- 
duced in Subsection 14. II The natural projection Q — > M x M is a principal SO(2)- 
bundle. As before, a curve y: [0, t] —> Q describes a rolling motion if there is no 
slipping, that is, if A(t)(x(t)) = x(t) and there is no twisting (see [2]) 

A(t) (vector field parallel along x(t)) — (vector field parallel along x(t)) . 

Let us now give expressions of the rolling distribution in local coordinates about 
a point (x,x;A) g Q. Let us consider local orthonormal frames ei,e2 forM and <?i,e2 
forM. They define their structure constants c\,C2 G C°°(M) and c\,C2 E C°°(M) by 
the equations [ci,C2] = c i e i +C2<?2 onM and [<?i,<?2] = c\e\ +C2e2 onM. 

Since Q is a circle bundle over M x M, in the natural trivialization, there a well 
defined angular direction and we can identify the isometry A with an angle 9. 
With these notations, the rolling distribution Z)r is spanned by the vector fields 

/ \ d 

X\ = e\ +cos 9ei + sin 9e2 + ( — c\ +c\ cos 9 + C2sin0) — -, 

o9 

, x B 

X2 = <?2 — sin0<*i +cos 9 £2 + I — C2 — c\ sin 9 + C2 cos 9 I ^r— . 

o9 

The main controllability theorem for the system of two Riemannian surfaces 
rolling, as presented in ||3] Chapter 24], is the following. 
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Theorem 4 (Agrachev & Sachkov 1999). Let G = G DR (q) be the orbit of the 
rolling distribution starting at q E Q and let K and k be the Gaussian curvatures of 
M and M respectively. Then: 

1. The orbit G is a imbedded connected submanifold of Q of dimension 2 or 5. 
More precisely, one has that if (k — k) | g is identically zero, then dim G = 2; 
and if (k— fc)\fi- is not identically zero, then dim G = 5. 

2. There is an injective correspondence between isometries 1: M — > M and two 
dimensional orbits of the rolling system. In particular, if the manifolds M and 
M are isometric, then the rolling model is not completely controllable. 

3. If M and M are complete and simply connected, then the correspondence be- 
tween isometries I : M — > M and two dimensional orbits of the rolling system is 
bijective. In particular, the rolling model is completely controllable if and only 
if the manifolds M and M are not isometric. 



5 Re-discovery of the higher dimensional case and interpolation 

Here we briefly review the way the higher dimensional problem of rolling manifolds 
presented to the control theory community and we explain how this was employed 
in geometric interpolation theory. 



5.1 Sharpe's definition 

Here we present the definition of rolling maps found in the Appendix B of R. W. 
Sharpe's book [48 1 with some minor modifications. 

Definition 4. Let M,M be n-dimensional submanifolds of IR" +V . Then, a differen- 
tiable map g : [0, f] — > Isom(R" +v ) satisfying the following conditions 

• There is a piece wise smooth curve x : [0, t] —> M, such that 

- g(t)x(t)€M, 

- T g {t) x {t) (g(t)M) = T gm) M. 

• Furthermore, the curve x(t) := g(t)x(t) satisfies the following conditions 

- No-slip: g(t)g(t)- 1 x(t)=0. 

- No-twist, tangential part: d{g{t)g{t)- l )T m M C T (g{t)g{t)- l M) L . 

- No-twist, normal part: d(g(t)g(t )- 1 )7> (f) M- L C T Q (g(t)g(t y l M). 

for any t E [0, t] is called a rolling map of M on M without slipping or twisting. 

We do not know whether Sharpe was aware of the existence of the paper B31 
at the time of the publication of his book, but his deduction of the "correct" defi- 
nition rolling maps follows the same structure as Nomizu's. Nevertheless, Sharpe 
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does obtain plenty of extra information. For example he shows that in the imbed- 
ded rolling system there is a deep relation with the Levi-Civita connections of the 
manifolds and the normal connections to the imbeddings. Besides this, he is able to 
prove precisely that rolling is transitive, that is 

Theorem 5 (Sharpe 1997). Let M ,Mi,M 2 C W +v be three n-dimensional sub- 
manifolds, such that they are tangent to each other at a common point p G Mq PI 
Mi (~lM 2 . Let y: [0, t] — > M be given such that y(0) = p. Assume that M\ rolls on 
Mo along the curve y with rolling map g\, and similarly let Mi roll on M\ along 
the curve J — giJ, with rolling map gi. Then M 2 rolls on Mq along the curve y, with 
rolling map g2gi and with image curve f = g2gl7 = g2j- 



5.2 Applications to geometric interpolation 

An interesting application of the rolling system has been in interpolation. The article 
where this idea appeared for the first time is [1331 for the case of the two dimensional 
sphere. Afterward it was extended successfully to arbitrary dimensional spheres, 
Grassmanians and to the special orthogonal groups in ||3"T| . This last application 
was employed in ||32l to study the motion planning of a rotating satellite. Later on 
in [30 1 the idea was also shown to work on Stiefel manifolds. 

The setting of the interpolation problem seems quite innocent. Let xo,X\ , . . . ,xn € 
M be measurements at times = to < h < • • • < tjv = T, and consider given initial 
and final velocities v e T XQ M and w £ T XN M. The interpolation problem consists in 
finding a C 2 curve y: [0, f] — > M satisfying 

y(/i) = Xi, y(0) = v, y(-r) = w, (Interp) 

and y minimizes the functional 




(Energy) 



Let £2 = {y e C 2 | y satisfies ( IInterpI )}. Then 
Theorem 6 (Crouch & Silva Leite 1991). Ify e Q minimizes (IEnergyI ). then 

D 3 (D \ 

^r+R[ Jt r,r)r=o, 

on each [f,-, f/+i], where R is the curvature tensor ofM. 

The curves satisfying the differential equation in Theorem [6] are called geomet- 
ric cubic splines, and they are in general quite hard to find. Nevertheless, in the 
cases described above, the authors were able to find a surprising relation between 
the rolling dynamics and geometric interpolation. For simplicity of exposition, we 
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only present the relevant results for the case of the n dimensional sphere S". A first 
observation that takes place is the following. 

Theorem 7 (Jupp & Kent 1987, Huper & Silva Leite 2007). Let R T (t) be the 
rotational part of the rolling map in Definition^ with rolling curve y: [0, t] — > S". 
For all t G [0, z] and all j G N, 

* T M^r(0 = rL +1) (0, 

where Ydev is the development of y, see Subsection \2.3\ 

A consequence of the above is the following application to interpolation in S". 

Corollary 1 (Jupp & Kent 1987, Huper & Silva Leite 2007). If the development 
1 1 y 7dev(0 is an Euclidean cubic spline, then 1 1-> 7(f) is a geometric cubic spline 
on S" if and only if it is a re-parameterized geodesic. 



6 Nowadays: The coordinate-free approach 

The intrinsic definition of the rolling model in higher dimensions was presented for 
the first time in 1 22, 26 1. It is clearly motivated by the definition given by Agrachev 
and Sachkov in iff]. 

Let (M,g) and (M,g) be two oriented Riemannian manifolds of dimension n. 
The state space of the rolling problem is the manifold 

Q = Q(M,M) = {A : T\ X M — > T\$M \x£M, x G M, Alinear isometry, det(A) >0}. 

An absolutely continuous curve q(t) = (y(f), y(t), A(f)) in Q is a rolling curve if 
A{t)X(t) is parallel along y(t) for every vector field X(t) that is parallel along 7(f) 
(no twist condition) and if A (f) 7(f) = f(t) (no slip condition). 

A counting argument shows that Q has dimension \n(n + 3). Over this manifold 
there is an n-dimensional distribution Dr, called the rolling distribution, such that 
the rolling curves in Q are exactly the integral curves of Dr. Let us describe this 
distribution briefly as given in 11221 . For a configuration q = (x,x;A) G Q, and an 
initial velocity X G T\ X M, we define the rolling lift _Sf R (X)| 9 G T\ q Q as 

^ R (X)\ q = j t \ (P' (y)°AoP?( 7 )), (7) 

where 7, y are any smooth curves in M, M, respectively, such that 7(0) = X and 
7(0) = AX, and (7) (resp. P*(y)) denotes the parallel transport along 7from y{a) 
to y(b) (resp. along 7 from y(a) to f(b)). 

Definition 5. (cf. Il22l ). The rolling distribution Dr on Q is the n-dimensional 
smooth distribution defined, for q = (x,x;A) G Q, by D R | 9 = ^f^(T\ x M)\ q . 
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An interpretation of the rolling lift Jz?r(X)| 9 of X e T\ X M at q = (x,x;A) is as 
follows. Let 7 be a curve in M such that 7(0) = x and y(0) = X then, by the general 
theory of ordinary differential equations, for short times there is a rolling curve q(t) 
of M on M satisfying q(0) = q. The rolling lift is precisely q(0). 



6.1 The controllability problem 

The orbit (q) of the rolling system described above passing through q 6 Q con- 
sists of all the states q that can be connected to q via a rolling curve. The (complete) 
controllability problem asks for conditions on the geometry of M and M such that 
ffo K {q) = Q. One way of addressing this problem is via Sussmann's orbit theorem, 
that is, by showing that all the Lie brackets of the vector fields steering the dynamics 
have to span the tangent bundle of the state space. For the rolling model, this Lie 
brackets are expressed in terms of the curvature tensors R and R associated to the 
Riemannian metrics g on M and g on M respectively, together with the covariant 
derivatives of R and R. It seems therefore impossible to solve for general dimen- 
sion n the controllability issue on the sole knowledge of the Lie algebraic structure 
of Dr, except for low dimensions. Indeed, in the case for instance where (M,g) 
is the n-dimensional Euclidean space, it would amount to determine Hol(V 4 '), the 
holonomy group of the Levi-Civita connection V 8 associated to g, with the only 
knowledge of its curvature tensor and its covariant derivatives. Instead, the latter 
issue can be successfully addressed by resorting on group theoretic and algebraic 
arguments, see [22|. For specific examples, using extra knowledge of the problem 
at hand, see E6ll36ll54l. 

In general, one can define a notion of curvature especially adapted to the rolling 
model, see l23l . For q — (x,Jc;A) e Q, the rolling curvature is the linear map 

2 

RoL,: /\T\ X M^ T*\ X M®T\ X M\ RoL,(X AY) := AR(X,Y) ~R(X,Y)A. 

This map permits to give a first sufficient condition for the rolling model to be 
controllable, see Il22ll28l . 

Theorem 8 (Chitour & Kokkonen 2011, Grong 2012). If the rolling curvature is 
an isomorphism for every q € Q, then the rolling problem is completely controllable. 

The above condition is very hard to deduce directly from conditions on the ge- 
ometry of M and M. It is therefore necessary to reduce the problem to a simpler one. 
One possible way to do this is to give some extra structure to the manifold M. In this 
vein, it was possible to give controllability conditions "without Lie brackets" for the 
case in which (M,g) = (F",g") is the space form of constant sectional curvature c, 
see Il37ll47l . To state these, let us first introduce some terminology. 

Definition 6. Consider the vector bundle Utakdm '■ TM ©I->M. The rolling con- 
nection V is the vector bundle connection on Utm®r defined by 
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VJ(y,a) = (v x Y + s(x)X,X(s)-cg(Y\ x ,X)), (8) 

for every x G M, X G T\ X M, (Y,s) G VF(M) x C°°(M); where we have canonically 
identified the space of smooth sections r(7iTM®M.) of 7Ztm®wl with VF(M) x C°°(M). 

When c^O, the connection V 6 is a metric connection with respect to the fiber 
inner product h c on TM R defined by 

h e ((X,r),(Y,s))=g(X ) Y)+c- 1 rs ) 

where X,Y E T\ X M, r,s G R. Its holonomy group is denoted by Jrf? c (M). In this 
language, we have the following result, see 11231 . 

Theorem 9 (Chitour & Kokkonen 2012). Let (M, g) be a complete, oriented and 
simply connected Riemannian manifold. The rolling problem of M rolling against 
F" is completely controllable if and only if 

( SO(n+l), c>0; 
je c (M) = \ SE(n), c = 0; 

[S0 (n,l), c<0. 

Here the Lie group SOo(n, 1) represents the identity component of the group 0(n, 1) 
of linear transformations that preserve the quadratic form F n \ (x\, . . . ,x n+ \) = x\ + 

2 2 
<~ x n ~ x n+V 

Wanting to fully understand these cases, it is important to remark some structure 
theorems encoded in Theorem [9] Observe that up to rescaling, it is sufficient to 
study when c = 0, 1 and —1. In the Euclidean situation, i.e. c = 0, the condition 
Jif°(M) = SE(n) is equivalent to that M has full Riemannian holonomy. In the case 
c = 1, if the action of Jf? 1 (M) on the unit sphere is not transitive, then (M,g) is the 
unit sphere. As a consequence, it holds that, for n > 16 and even, the rolling system 
Q = Q(M,S"~ l ) is completely controllable if and only if (M,g) is not isometric to 
the unit sphere. Both theses cases were analyzed in [23 1, and the remaining cases 
are currently under investigation. The hyperbolic case presented a more difficult 
challenge, see fl9l . 

Theorem 10 (Chitour, Godoy & Kokkonen 2012). Let (M, g) be a complete, ori- 
ented and simply connected Riemannian n-manifold rolling onto the space form 
(H",g" j) of curvature —I. Then the associated rolling problem is completely con- 
trollable if and only if (M,g) is not isometric to a warped product of the form 

(WP1) (RxMi,ds 2 ®e*gi), or 

(WP2) (H* x Mi , gij © cosh (^ rf) g\ ), where \<k<n and for each x G , d{x) 
is the distance between x and an arbitrary fixed point xq G H . 

In both situations, (M\,gi) is some complete simply connected Riemannian mani- 
fold. As usual, the term ds 2 represents the usual Riemannian metric on R. 
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6.2 Symmetries of the rolling problem 

The idea developed in Subsection l6.1l of setting M to be a space form has a beautiful 
geometric consequence on the bundle structure of the natural projection Kq^m '■ Q = 
g(M,F") — > M. Let us explain what this is. 

In general, it is not clear if there is a G-principal bundle structure on Q making 
Dr a G-principal bundle connection. This is indeed the case if the manifolds are of 
dimension 2, in which case the projection Q — >• M x M is a principal SO(2) bundle 
with Dr as its connection. Nevertheless, for higher dimensions the projection Q —> 
M xM does not satisfy the above. The main reason is that SO(n) is abelian only for 
n = 2, thus we need to make the problem simpler. 

For c^O, let G c (n) be the identity component of the Lie group of linear 
maps that leave invariant the bilinear form (x,y)" := Yli=\ x iyi + c X n +iy n +\, for 
x = (x\,... ,x n+ i),y — (yi,... ,y„+i) S K n+1 . Observe that G\ (n) = SO(« + 1) and 
G_i(n) = SOo(«, 1). For c — 0, we set Go(«) = SE(«). Recall that, with this nota- 
tion, the identity component of the isometry group of (F",g") is equal to G c (n) for 
all c € K (cf. E3). 

The fundamental result concerning rolling against a space form lies in the fact 
that there is a G e (n)-principal bundle structure for the state space compatible with 
the distribution Dr, i.e. Dr is a G 6 (n)-principal bundle connection, see [23|. The 
precise result follows. 

Theorem 11 (Chitour & Kokkonen 2012). Let Q = Q(M, F") be the state space of 
rolling M against the space form F". Then we have: 

(i) The projection Kq^m '■ Q ~^ M is a principal G c (n)-bundle with a left action /J. : 
G c (n) x Q — > Q defined for every q = (x,x;A) by 

H(($,C),q) = (x,Cx+y;CoA),ifc = 0,n(B,q) = (x,Bx;B oA),ifc ^ 0. 

Moreover, the action preserves the distribution Dr, i.e., for any q G Q and 
B G G c {n), (p B )*D R \ q = D R \^ {Bq) , where pL B :Q^Q; q^ ll{B,q). 

(ii) For any given q — (x,x;A) 6 Q, there is a unique subgroup ,JP q c of G c (n), called 
the holonomy group of Dr at q G Q, such that 

M(^ C x {?}) = 0Di(q) n 7t Q l M (x). 

Also, if q' = (x,j?;A') G Q is in the same %Q M~fiber as q, then ,^ q c and $? c , 
are conjugate in G c (n) and all conjugacy classes of M'^ in G c (n) are of the 
form .y^f"). 

J q' 

A natural question to ask is whether a converse of the theorem above holds, in 
other words, does the existence of a G-principal bundle structure on Q such that Dr 
is a connection imply that M must have constant sectional curvature? The answer is 
generically yes, but we need to introduce some more terminology. 

Recall that in Subsection l2.2l we defined the Lie algebra of symmetries Sym(D) 
of a distribution Dona manifold M as the set of vector fields X G VF(M) that 
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satisfy [X,D] CD. For the case of the rolling distribution, we will focus our attention 
in the symmetries of the rolling distribution that are annihilated by the projection 
Kqm '■ Q—t M, that is, in the Lie algebra 

Sym (D R ) := {S e Sym(D R ) | (tz qm ),S = 0}. 

With this at hand, the mentioned converse takes the following form, see [20|. 

Theorem 12 (Chitour, Godoy & Kokkonen 2012). If there is an open dense set 
O C Q such that R\ x : /\ 2 T\ X M -> /\ 2 T\ X M is invertible on Kq M {0) and Rol is 
invertible on O, then, up to an isomorphism of Lie-algebras, 

Sym (D R )=Iso(M,|) 

and therefore all the elements o/Symo(Z) R ) are induced by Killing fields of (M,g). 

In particular, under the above assumptions, if there is a principal bundle struc- 
ture on %q,m : Q—> M that renders Z) R to a principal bundle connection, then (M,g) 
is a space of constant curvature. 



6.3 Generalizations and perspectives 

Two natural questions to ask concern the extension of the rolling system to the 
situation in which the manifolds involved have different dimension and to to extend 
the classification result in Subsection |43] to other cases. For the first question, one 
needs to consider curves of isometric injections instead of isometries. This change 
introduces many difficulties in understanding the controllability problem, and in 
fact many tools that work well in the classical situation can not be generalized. The 
second question has a satisfactory answer for the three dimensional case, see [22 1 . 
There it is shown that the orbits can have dimensions 3, 6, 7, 8 and 9. 

A question that has been in our minds for a while is to actually compare the man- 
ifolds via the rolling model. This idea of comparison is naively evident in the rolling 
curvature tensor: one is actually subtracting the Riemannian curvatures of the man- 
ifolds. In fact, rolling should provide a framework for the isometric characterization 
of manifolds by using curvature tensor spectrum information (as in Osserman-type 
conditions for instance, cf. [25]). 

Finally, we have noticed that the problem of rolling manifolds can be generalized 
far beyond than allowing arbitrary connections, as in 1128] Section 7], or to pseudo- 
Riemannian manifolds, as in | : 4 I |. This extension consists of rolling so-called Car- 
tan geometries, see [48 1, and it includes as particular cases both of the situations 
mentioned above, together with the problem of rolling manifolds of different di- 
mensions, see [21 J. The main idea behind this is that Cartan geometries are the 
most general framework for a notion of development to exists, which underlies the 
very definition of the rolling dynamics. So far this generalized model has resisted a 
thorough study of controllability. 
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